NPS  ARCHIVE 
1958 
JOHNSON,  F. 


'  '  ;!!■'      ''       !l 


AN  IBM-701  COMPUTER  ANALYSIS 
OF  A  RESERVOIR  BY  MEANS 
OF  THE  MUSKAT  EQUATION 

FRANCIS  C.  JOHNSON 


DUDLEY  KNOX  LIBRARY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA   93943-5101 


An  IBM- 701  Computer  Analysis 

of  a  Reservoir  by  Means  of  the 

Muskat  Equation 


by 
Francis  C.   Johnson 


Petroleum  Engineering  299 
University  of  California 
Berkeley,    California 
Spring  Term,   1958 


DUDLEY  KNOX  LIBRARY 
NAVAL  POSTGRADUATE  SCHOOl 
MONTRREY  CA   93943-5101 


ABSTRACT 


A  method  is  presented  for  the  expeditious  analysis  of  performance 
of  a  solution  gas  drive  reservoir.     This  rapid  analysis  is  accomplished  by 
the  repeated  solution  of  the  Muskat  equation  on  the  IBM-701  computer.     In 
each  solution,  one  or  more  reservoir  parameters  are  varied  to  determine 
the  effect  on  the  problem  results.     A  small  California  reservoir  is  used  to 
demonstrate  the  procedure.     Graphs,   illustrating  the  results  for  each 
parameter  range,   are  presented. 

For  the  sample  reservoir  only  production,  depth-pressure,   and 
electric  log  data  were  available.     Average  values  of  PVT,   viscosity,   and 
permeability  were  developed  from  published  correlation  in  California  Oils, 
and  entered  into  the  Muskat  prediction  computer  program.  *    Fifty  predictions 
were  run,   and,  of  these,  twenty-three  solutions  are  graphically  illustrated. 
From  the  graphs,   a  set  of  parameters  was  selected  which  best  define  the 
sample  reservoir. 

The  results  of  this  analysis  indicate  that  more  than  one  prediction 
is  necessary  to  establish  the  reservoir  parameters.     By  continuously  ob- 
serving the  predictions,  the  reservoir  engineer  is  able  to  decide  when  a 
reasonable  approximation  of  the  reservoir's  parameters  has  been  made. 
He  must  know,  from  experience  and  his  knowledge  of  reservoirs,   which 
set  of  parameters  best  defines  the  reservoir. 
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INTRODUCTION 

For  the  last  decade,    Muskat's  equation'  '  for  the  depletion  history 
of  a  solution  gas -drive  reservoir  has  been  a  reliable  tool  of  the  reservoir 
engineer.     It  undoubtedly  would  be  used  more  extensively,  but  for  the  vast    . 
amount  of  time  and  labor  required  to  complete  one  Muskat  reservoir  pre- 
diction.    (An  estimated  sixteen  man  hours  per  prediction).     Often,   some  of 
the  data  which  the  reservoir  engineer  must  use  are  doubtful.     It  is  desirable, 
therefore,   that  the  Muskat  prediction  be  run  a  number  of  times,   using  an 
array  of  values  for  the  doubtful  data.     However,  the  man  hours  required  for 
such  a  procedure  would  soon  become  prohibitive.     The  obvious  solution  to 
this  dilemma  is  the  high  speed  digital  computer.     Muskat's  material  balance 
equation,   expressed  in  differential  form,  is  ideally  suited  to  digital  computer 
solution. 

The  purpose  of  this  paper  is  to  illustrate  how  a  series  of  Muskat  pre- 
dictions, obtained  in  a  relatively  short  time,  provided  the  reservoir  engineer 
with  a  production  analysis  of  his  reservoir. 

A  small  California  reservoir  was  used  to  demonstrate  this  procedure. 
Data  for  the  reservoir  were  provided  by  the  Standard  Oil  Company  of  California. 
These  data  included  electric  logs,  production  records,  and  a  limited  amount 
of  depth  pressure  data.    No  PVT  data  or  core  data  were  available.     Porosity 
and  water  saturation  data  for  a  similar  zone,   in  a  nearby  reservoir,  were 
used  to  make  initial  estimates. 
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THEORY 

The  IBM- 701  computer  solution  to  the  Muskat  Prediction,   once  pro- 
grammed,  may  be  run,   with  slight  alteration  to  the  variables,    until  the  pre- 
diction data  and  the  production  data  coincide.     When  this  occurs,   it  may  be 
assumed  that  the  prediction,   extended  beyond  the  current  produced  life  of 
the  reservoir,   will  be  an  accurate  estimate  of  the  reservoir's  future,   pro- 
vided the  operating  procedure  of  the  reservoir  is  unchanged.     The  computer 
makes  a  single  prediction  in  about  three  minutes.     However,   if  three  pre- 
dictions are  run  in  tandem,   with  the  same  imput  quantity  varied  in  each  (e.  g„ 
a  different  size  gas  cap  in  each  run),  the  three  predictions  will  be  completed 
in  less  than  seven  minutes. 

The  IBM- 701  program  herein  utilized  was  developed  by  J.   E.   Warren 
and  T.   D.    Mueller2  of  the  California  Research  Corporation  at  La  Habra, 
California.     The  program  includes  the  equations  listed  below: 
For  production  above  the  bubble  point  pressure: 

1.     NP=   7758*S0|{^--fe[l  +  M,(l--p-)]] 

RSNP 
2*    GP  =     1000 

For  production  below  bubble  point  pressure: 

d_So_   SpQ  +  U  (I-S0-Sw)  +  SqY  k9/k0+M   BoSoi*QtU(l-Sw-BoSoi*)-Soi*fp-0 

' dp "  >+r%9)(Vko) 

^  '  d  Rs 

with     Q=  5.615  XB0      d~p~ 

u      r     dp 


3. 


X=/B0 


Y  = 


x      dBp 
XB02  dp 


4.  R  =  Rs+  5.615  X  k^Q 

5.  Np  .7758*(g--fl) 

6.  Gp  =   43.56  ♦  jtslllji  "    57IIJ,  +   ^  Sgi  -  /  Sg] 

.„[„.[,    S0i*   _    \     vf.    So;*        \  .   RsiSpi        RsSpj* 
+  ML/.II-   g--Sw)-yU-   b3     S«>+ 5.615  B0j      5.615  Bo 


where: 

N     *  Cumulative  oil  production  in  bbl  of  STO/acre  ft. 

G     *  Cumulative  gas  production  in  MCF/acre  ft, 

R     *  Produced  gas -oil  ratio  in  SCF/STbbl. 

All  other  symbols  are  defined  in  the  appendix,   page  i. 

The  computer  program  is  designed  to  solve  these  equations,   and 
to  record  their  solution  as  illustrated  by  Figure  1. 

To  adapt  the  computer  program  to  a  particular  reservoir,  the  PVT 
data,   viscosity  data,   and  permeability  data  for  that  reservoir  are  introduced 
into  the  program.     To  consider  a  new  reservoir,  one  simply  substitutes  new 
data. 
RESERVOIR  HISTORY 

From  the  electric  logs  provided,   a  sand  count  was  made  of  the  pro- 
ducing sand,   and  a  contour  map  was  developed,     (see  Figure  2).     It  was 
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apparent,  from  the  production  data  and  from  numerous  cross  sections  of 
the  producing  zone  and  the  two  zones  above  it.   that  the  reservoir  provided 
was  actually  three  reservoirs,   separated  from  one  another  by  fault  zones. 
It  was  decided  that  only  one  of  these  reservoirs.   Reservoir  B,   (illustrated 
in  Figure  2)  would  be  used  to  demonstrate  the  IBM- 701  prediction  procedure. 

The  contour  map  of  Reservoir  B  shows  that  the  reservoir  contains 
nine  wells.     Of  the  nine  wells  drilled,   only  two  are  currently  producing. 
Well  #5  was  the  first  well  brought  in  (November  1949)  with  an  initial  pro- 
duction of  400  bbls/D.     It  produced  until  December  1950  (a  total  of  82,  000 
bbls),   when  it  was  shut  in  with  a  maximum  gas-oil  ratio  of  7500.     Well  #18 
was  completed  in  January  1950,   initial  production  9  bbls/D.     It  produced 
until  December,   1951  when  its  maximum  gas-oil  ratio  was  6700o     It  produced 
a  total  of  2,  000  bbls  of  oil.     Well  #1  was  drilled  in  June  1950,   and  came  in  a 
gas  well.     It  was  immediately  shut-in.     Well  #8  was  completed  in  November 
1950;  its  initial  production  was  380  bbls/D.     It  is  currently  producing  31  bbls/D. 
Its  gas -oil  ratio  has  climbed  from  280  to  4,  000,     It  has  produced  a  total  of 
185,  000  bbls.   of  oil.     Well  #9  was  drilled  in  January  1953.     It  came  in  a  gas 
well  and  was  shut  in  after  producing  3,  600  MCF  of  gas  and  33  bbls.   of  oil. 
Well  #12  was  completed  in  January,   1956  with  an  initial  production  of  49  bbls/D. 
It  has  produced  22,  800  bbls.   to  date  and  is  currently  producing  51  bbls/D.     This 
well  initially  produced  118  bbls/D  of  water,    and  is  now  producing  272  bbls/D. 
It  is  the  only  well  in  the  field  that  has  produced  a  sizable  amount  of  water. 
Wells  #2,   #4,   and  #14  were  dry  holes.     As  illustrated  in  Figure  2,   the  entire 
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northern  part  of  Reservoir  B,   which  includes  wells  #5,    #18,    #1,    #2,   and  #9 
R/D,   is  a  gas  cap.     There  is  also  evidence  of  water  influx  in  the  eastern 
edge  of  the  reservoir  near  well  #12. 

The  total  oil  produced  to  date  from  Reservoir  B  (July  195  7)  is  507,  356 
bbls.   of  stock  tank  oil.     This  oil  has  an  average  gravity  of  35°  API.     673,  457 
MCF  of  gas  has  been  produced. 
PREDICTION  PROCEDURE 

Production  data  provided  an  average  API  gravity  of  35°  for  stock 
tank  oil,   an  average  gas  gravity  of  0.  73  for  the  reservoir  gas  at  standard 

stock  tank  conditions,   and  an  average  reservoir  temperature  of  135°  F. 

3  4 

Standing's     correlation  charts  were  used  to  obtain  PVT  data,   and  Beal's 

correlations  were  used  for  the  viscosity  data,     PVT  data  and  viscosity  data, 

as  a  function  of  pressure,   are  included  as  Graph  A. 

From  the  sand  count  of  Reservoir  B,   an  isopachous  map  was  constructed 

(Figure  3).     A  volumetric  estimate  of  sand  in  place  was  computed  to  be  2,  520 

acre  feet.     Porosity  and  connate  water  saturation  were  estimated  to  be  0.  35 

and  0.  45  respectively.     These  were  the  average  values  of  twenty-two  porosities 

and  water  saturations  taken  from  sand  in  an  adjacent  reservoir.     The  adjacent 

sand  was  considered  to  be  similar.     Using  the  equation: 
N    -      (7756)  (acre  ft.)  (<j>  )  (1  -  S  w)      , 

Boi 
an  initial  estimate  of  oil  in  place  was  computed  to  be  2.  88  million  stock  tank 

barrels.     (Sample  calculation,   page  ii  in  appendix  illustrates  this  procedure). 
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This  estimate  was  substantiated  by  the  Schilthuis  Material  Balance  Equation 
when  an  allowance  was  made  for  gas  production  from  the  gas  cap.     A  second 
calculation  of  N  was  made  using  new  values  of  porosity  and  connate  water 
saturation,   <t>  =   0.  25  and  Sw  =0.  30.     These  estimates  were  taken  as  average 
for  California  sandstone.     The  initial  estimates  are  on  the  extreme  end  of 
the  range  for  California  sandstones.     The  new  values  of  porosity  and  water 
saturation  give  a  value  of  N  equal  to  3. 14  million  barrels  of  stock  tank  oil. 
The  change  in  N  is  comparatively  slight. 

A  curve  of  oil  saturation  vs.   relative  permeability,   (Graph  B),    was 
developed  by  solving  the  equations  listed  below: 

7.     So.{l-Sw)(l-^)g 


8.     k9/k0   =  (R-Rs)M9/Vo7b0 


Static  bottom  hole  pressure,   cumulative  oil  and  gas  production,  and 
instantaneous  gas -oil  ratios  were  plotted  vs.   time  in  years,     (see  Graphs  C 
&  D).     From  the  correlations  established,   R  and    AN  (bbls.   of  STO)  were 
taken  as  functions  of  pressure.     Dividing    A  N  by  N  gave     A  N/N  as  a  function 
of  pressure.     Since  two  N's  were  computed,  two  cases,   Case  I,   where  <t>  ■ 
0.  35  and  Sw  *  0.  45  and  Case  II,   where  <t>  =  0.  25  and  S     s  0.  30  were  developed. 
PVT  and  viscosity  data,   as  functions  of  pressure,   were  taken  from  Graph  A. 
Equations  (7)  and  (8)  were  solved  as  functions  of  pressure  in  decrements  of 
100  psi  from  1500  to  600  psia.     From  the  correlation  of  relative  permeability 
and  oil  saturation  for  the  same  decrement  of  pressure,   two  curves  of  K^/Kq 
vs.  S0  for  Case  I  and  Case  II.     Both  curves  were  extended  into  the  region  of 
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higher  gas  saturation  by  continuing  them  parallel  to  the  curve  developed 
by  Arps     for  the  average  value  of  relative  permeability  in  sandstones.     Once 
the  PVT,   viscosity,  and  relative  permeability  data  were  compiled,   they  were 
included  in  the  computer  program.     Experience  with  the  computer  program 
indicates  that  most  functional  data,    used  by  the  computer,   must  be  mathema- 
tically smoothed.     The  smoothing  of  the  PVT  and  viscosity  data  was  accom- 
plished by  third  order  polynomial  equations.     A  seventh  order  equation  is 
usually  required  for  K_/K0  functions.     However,   in  this  particular  case, 
the  curve  from  Arp's  average  was  sufficiently  smooth  to  enter  the  program 
directly. 
PREDICTION  RESULTS 

It  was  anticipated  that  the  predicted  results  could  be  brought  to 
superimpose  upon  the  observed  production  data.     This  was  to  be  accomplished 
by  a  system  of  altering  the  problem  variable  until  the  desired  results  were 
reached. 

For  the  first  prediction,  it  was  decided  to  investigate  the  Case  I 
parameters.     The  following  initial  conditions  were  set: 

Pj    =  1500  psia,  bubble  point  pressure  =  1500  psia 

Sw    =   0.45,   <f>  =  0.35,   m  «   0.1,   S0i*  «   0.15 
The  predicted  results  are  illustrated  in  Figures  4.  0  and  4.1.  * 


♦Predicted  data  was  illustrated  in  all  figures  by  a  solid  curve.     Observed 
data  was  shown  as  circled  points. 
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Figures  5.  0  and  5.1  show  that  lowering  the  bubble  point  causes  the 
predicted  data  to  begin  to  approach  the  observed  data.     The  best  approxima- 
tion of  %     A  N/N  is  for  a  BPP»  900  psia,  but  the  declining  slope  after  P* 
900  is  too  gentle. 

After  this  prediction  was  made,   the  effect  of  varying  porosity,   water 
saturation,   viscosity,   gas  cap  and  initial  pressure  was  investigated.     None 
of  these  variations  produced  a  unique  solution.     Table  1.  0  lists  the  parameters 
investigated. 
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Very  little  agreement  is  shown  between  observed  and  predicted 


data. 


To  investigate  the  effect  of  bubble  point  on  the  prediction  a  second 
run  was  made.     This  was  a  tandem  run  with  a  different  bubble  point  used 
for  each  prediction.     The  following  initial  conditions  were  set: 

Pi  ■  1500  psia,BPP  *  1500,   1200,   900  psia,   Sw  »  0,  45,   <t>  ■   0.35, 

m  =   0.1,   Soi  *  0.15 
The  predicted  results  are  illustrated  in  Figures  5.  0  and  5.1. 
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TABLE  1.0 


CASE  I 


Parameters  Not  Illustrated  on  Graphs 


Initial 

1500 

Pressure 

1500 

1500 

1500 

1500 

1500 

1400 
1300 

1500 

1300 

Bubble 

1100 

Point 

1500 

900 

900 

900 

900 

1500 

900 

1000 

Pressure 

900 

0.30 

Porosity 

0.  35 

0.35 

0.35 

0.35 
0.40 

0.35 

0.  35 

0.  35 

0.  35 

Connate 

0.40 

Water 

0.45 

0.45 

0.45 

0.45 

0.45 

0.45 

0.45 

0.45 

Saturation 

0.50 

Gas  Cap  Vol. 

0.1 

0.1 

Oil  Vol. 

0.5 

0.5 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

(m) 

0.9 

0.9 

M 

0.  8  M 

Viscosity 

1.0 

1.0 

1.0  M 

1.0 

1.0 

1.0 

1.0 

1.0 

Ratio 

1.  5  M 

0.  20 

S   •* 
°oi 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 
0.10 

0.15 

28.6 

28.1 

30.4 

26.4 

28.  6 

28.0 

28.5 

%    AN/N 

33.1 

33.  6 

28.1 

28.1 

28.1 

27.  8 

28.1 

28.1 

36.  2 

37.  2 

24.  2 

30.1 

26.9 

28.3 

25.4 

The  Case  II  parameters  were  next  investigated.     For  the  third  pre- 
diction, the  following  initial  conditions  were  set: 

PA    «  1500  psia,  BPP    ■  1500  psia,   Sw    *  0.30,   <j>    *  0.  25,   m    =   0.5, 
0.15 
The  predicted  results  are  illustrated  in  Figures  6.  0  &  6.1. 
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Here  the  GOR  curve  showed  slight  agreement.  No  agreement  was 
apparent  for  oil  production. 

The  fourth  prediction  was  used  to  investigate  the  effect  of  gas  cap 
size.     The  following  initial  conditions  were  set: 

Pi    -  1500.  BPP  -  1500  psia,   Sw    -   0.30,   <t>    *   0.  25,   m    =  0.1,0.5 

&  0.9,   Soi*    »   0.15 


The  predicted  results  are  illustrated  in  Figures  7.  0  and  7. 1. 
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The  fourth  prediction  indicated  that  the  gas  cap  size  does  have  a 
significant  effect  upon  gas  and  oil  produced.     It  suggests  that    m  *  0. 1  is 
probably  the  more  likely  size  for  this  reservoir.     There  was  no  agreement 
between  predicted  and  observed  data  for  oil  production,  but  predicted  and 
observed  gas -oil  ratio  compared  well. 

The  fifth  prediction  was  used  to  demonstrate  the  effect  of  bubble 
point  upon  production.     For  this  prediction  the  following  initial  conditions 
were  set: 

Pi    *  1500  psia,  BPP    «  1500,    1200,    900  psia,   Sw    »   0.  30,   <t>    *   0.  25 

m    =0.1,   Soi*    *   0.15 


The  predicted  results  are  illustrated  in  Figures  8.0  and  8.1- 
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The  fifth  prediction  showed  that  lowering  the  bubble  point  pressure 
caused  the  predicted  oil  curve  to  approach  the  curve  of  the  actual  oil  produced. 
At  a  bubble  point  of  900  psia,  the  predicted  oil  curve  superimposed  upon  the 
actual  produced  oil  curve.     However,   the  corresponding  GOR  curve  at  BPP  » 
900  psia  was  too  low.     This  may  be  explained  by  examining  the  history  of 
Reservoir  B.     The  excessively  high  gas  oil  ratios  which  occur  at  various 
times  throughout  the  reservoir's  history,  strongly  suggest  that  gas  was 
produced  from  the  gas  cap.     Since  the  Muskat  equation  does  not  account 
for  gas  production  directly  from  the  gas  cap,   the  computer  program  was 
not  designed  to  reflect  it.     It  is  concluded  that  the  difference  between  GOR 
computed  and  actual  GOR  is  due  to  the  gas  produced  directly  from  the  gas 
cap.     The  prediction  indicates  that  slightly  over  half  of  the  gas  produced 
from  the  initial  pressure  of  1500  psia  to  600  psia  was  produced  directly 
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from  the  gas  cap.     A  new  estimate  of  actual  GOR  curve  was  made  excluding 
the  gas  cap  production.     This  new  GOR  will  be  reflected  on  the  remaining 
curves.     The  fifth  prediction  suggests  that  a  bubble  point  of  approximately 
900  psia  is  most  likely  for  Reservoir  B.     Using  this  new  bubble  point  results 
in  a  new  value  of  N,  based  on  a  change  in  B0i-     The  change  is  slight,   however, 
and  not  sufficient  to  alter  the  Kg/K0  vs.  S0  curve. 

It  should  be  noted  on  figure  8.  0  that  the  recoverable  oil  at  abandonment 
pressure  of  200  psia  was  roughly  the  same  for  all  three  bubble  points. 

To  substantiate  the  earlier  decision  to  use  a  gas  cap  size    m  *  0.1, 
in  the  sixth  prediction,  a  variable  gas  cap  was  once  more  employed;  this 
time  with  the  bubble  point  pressure  at  900  psia.     The  following  initial 
conditions  were  set: 

Pi  -  1500  psia,  BPP«  900  psia,  Sw  *  0.30 

4>   «   0.  25,     m  -   0. 1,    0.  5,   &  0.  9,   Soi  =  0. 15 
The  predicted  results  are  illustrated  in  Figures  9.  0  and  9.1. 
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Once  again,   an    m  =  0.1    appeared  to  be  the  best  estimate  for  gas 
cap  size  in  Reservoir  B. 


To  investigate  the  influence  of  viscosity,  the  seventh  prediction  was 

made  with  the  following  initial  conditions: 

Pj    «  1500  psia,   BPP  *  900  psia,   Sw  «  0.  30,   <t>    *   0.  25, 

m    ■0.1.  Soi*  -  0. 15,  and  0.  8  M,  1.  OM  &  1.  5M  where  M  «  M*Im<1 

The  predicted  results  are  illustrated  in  Figures  10.  0  and  10. 1. 
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As  one  would  expect  a  higher  viscosity  results  in  less  recoverable 
oil  and  a  consequent  higher  GOR.     The  results  indicated  that  the  initial 
estimate  of  oil  viscosity  (  M    *  ^'iMf    x  1.  0  )  is  the  best  approximation  for 
the  actual  production  data. 


In  the  eighth  prediction,   the  influence  of  porosity  was  investigated. 
Changing  porosity  changed  the  value  of  N,  the  initial  oil  in  place.     It  did  not 
change  the  value  of  %     A  N/N  because  porosity  cancels  out  of  this  ratio. 
However,   in  the  calculation  of  oil  and  gas  produced  per  acre  foot,   porosity 
is  reflected.     Converting  the  acre  feet  production  to  predicted  reservoir 
production  permits  a  comparison  with  actual  production.  Graph  D. 

The  change  in  N  for  each  value  of  porosity  resulted  in  a  very  slight 
shift  of  Kg/KQ  along  the  SQ  axis  (Graph  B  and  Equation  7).     The  shift  was 
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so  slight  that  Arp's  average  curve  was  not  varied. 

The  eighth  prediction  initial  conditions  were  set  as  follows: 

Pt    -  1500  psia,  BPP    «   900  psia,   Sw    »   0.  30,     <t>    -   0.  20,    0.  25  &  0.  30. 

m    *   0.1,    Soi*    *   0.15 

The  predicted  results  are  illustrated  in  Figures  11.  0  and  11. 1. 
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This  prediction  is  compared  with  actual  production  data  illustrated  on 
Graph  D.     The  production  data  is  not  altered  by  dividing  it  by  N.     It  is  concluded 
that  for  Reservoir  B,   a  porosity  of  0.  25  is  the  only  porosity  which  will  give  a 
prediction  that  satisfies  actual  production  data.     This  condition  holds  so  long 
as  Kg/K0  curve  was  shifted  along  the  SQ  axis  for  each  change  in  Sw.     This  shift 
was  not  due  to  change  in  N,  but  due  entirely  to  the  change  in  Sw. 


The  ninth  prediction  was  made  to  investigate  the  effect  of  connate  water 
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saturation.     With  the  porosity  held  constant  at  0.  25,   the  change  in  Sw  caused 
a  slight  change  in  N.     To  illustrate  how  this  change  in  N  caused  a  shift  of 
Kff/K     curve  along  the  S0  axis,   a  sample  calculation  was  made.     N  is  expressed 
in  millions  of  stock  tank  barrels.     When  positioning  the  Kg/ Kq  curve  on  the 
S0  axis,   S0  was  obtained  by  solving  Equation  (7)  at  various  values  of  pressure. 
Example:    when  P  *   700  psia, 

S0  -  (1  -  Sw)    (1  -    ^-     ) 


N  Bi 

at    N  '   3. 12    when  <t>  =  0.  25  and    Sw   *   .  30 


SQ  =  (.7)    (1  -  l|1      )    (.865)     x   0.56 

fill 


if  N    =    2.  9    for  <j>  =   0.  25  and  Sw  *   3.5 


S0  *  (.  65)  (1  -    -41      )    (.  865)    =   0.  515  ,    a  shift  of  4.  5%. 
2.9 

However,   if  only  the  change  in  N  is  reflected, 
N   '   2.  9  for  <j>  =   0.  25  and  Sw  =   0.  30 

SQ  *  (.7)    (1  -    ^-gp)    (.865)    =   .555 

or  the  shift  of  the  Kff/KQ  curve  along  the  SQ  axis  for  a  change  in  N  of  220,  000 
bbls.   of  STO  is  only  0.  5%.     It  was  concluded  that  the  small  changes  in  N,   due 
to  changes  in  porosity  and  water  saturation,    were  not  sufficient  to  shift  the 
Kg/KQ  curve,  but  the  changes  in  Sw  did  move  the  curve  along  the  S0  axis  a 
distance  equal  to  that  change.     Consequently,   it  follows  that  Kg/K    vs.  S0 
curve,  taken  as  Arp's  average  for  sandstones,   is  a  fixed  parameter  except 
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when  the  value  of  Sw  is  varied. 

Initial  conditions  for  the  variable  connate  water  saturation  prediction 


were  set  as  follows: 

Pi  *   1500  psia,  BPP  =   900  psia,    Sw  -   0.  25,    0.  30,    and  0.  35,   <t> 

m  *   0.1,    S0i*  *   0.15. 
The  predicted  results  are  illustrated  in  Figures  12.0  and  12.1. 
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Figures  12.0  and  12.1  illustrate  the  effect  of  the  new  Kg/KQ  vs.  S0 
curves  on  the  prediction.     These  curves,   which  were  extended  by  Arp's 
average,   satisfied  the  observed  data  only  when  Sw  =   0.  30. 


The  tenth  prediction  was  made  to  investigate  the  effect  of  initial 
pressure.     Initial  conditions  were  set  as  follows: 

P.  »  1500,   1400,   1300  psia.  BPP  =  900  psia,   Sw  =   0.30,   6  =   0.25, 
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m  -   0.1,    Soi*  »   0.15. 
The  results  were  not  of  sufficient  difference  to  merit  plotting.     The  oil 
recovery  at  abandonment  pressure  of  200  psia  for  each  Pj  was: 
Pt  *  1500  psia.   %     AN/N    «   23.57 
1400  psia,    %      AN/N    -   23.49 
1300  psia,   %      AN/N    *   23.39 
It  was  concluded  that  an  initial  pressure  from  1500  to  1300  psia  will 
have  only  slight  effect  on  recoverable  oil. 

The  eleventh  prediction  illustrated  the  effect  of  changing  S0j*.     The 
following  initial  conditions  were  set: 

Pi  *  1500  psia,  BPP  =   900  psia,   Sw  =   0.  30,    <t>  »   0.  25 

m  =   0.1,   Soi*  =   0. 10,    0. 15,   and  0.  20. 

The  results  indicated  that,   within  the  average  limits  of  S0j*  for  a 
California  solution-gas  drive  reservoir,   S0i  variation  has  little  or  no  effect. 

It  was  concluded  that  the  following  parameters  describe  Reservoir  B: 

Bubble  point  pressure 900  psia 

Initial  reservoir  pressure 1300  -  1500  psia 

Connate  water  saturation 0.  30 

Porosity 0.  25 

Ratio  of  initial  reservoir  free  gas  volume  to  initial  reservoir  oil 
volume 0.1 
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Oil  saturation  in  the  gas  cap 0. 15 

PVT,    viscosity,   and  relative  permeability  -  -  -  as  described  in  Graphs 
A  and  B,    respectively 

Initial  oil  in  place  in  Reservoir  B  was  estimated  to  be  3.12  million 
barrels.     Recoverable  oil  was  estimated  to  be  730,  000  barrels  of  stock  tank 
oil  or  23.  57%  of  initial  oil  in  place. 
DISCUSSION  OF  RESULTS 

The  prediction  results*    have  illustrated  how  quickly  the  reservoir 
engineer  may  bring  predicted  data  to  superimpose  upon  actual  production 
data.     By  the  fifth  production,   it  was  reasonably  apparent  that  recoverable 
oil  was  about  23.  6%  of  oil  in  place.     The  remaining  predictions  simply  re- 
affirmed the  fifth  prediction  by  illustrating  that  the  Case  II  values,   selected 
for  the  variable  data,   gave  solutions  which  more  nearly  approximated  the 
actual  production  data.     Such  will  not  always  be  the  case.     However,   when 
the  reservoir  engineer  runs  the  Muskat  prediction  in  tandem,    with  three 
different  values  for  some  doubtful  input,   he  will  very  rapidly  pick  up  a 
trend  which  will  enable  him  to  close  in  on  the  parameters  which  specify  his 
reservoir. 

Reservoir  B  was  not  the  best  possible  reservoir  to  illustrate  this 
problem.     It  is,   however,   typical  of  the  situations  which  the  reservoir 
engineer  continuously  encounters. 

From  the  history  of  Reservoir  B,   it  was  obvious  that  gas  was  often 
*See  Table  2.  0. 
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TABLE  2.0 


Summary  of  Parameters  Used 


Prediction 
number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Initial 
Pressure 

1500 

1500 

1500 

1500 

1500 

1500 

1500 

1500 

1500 

1500 
1400 
1300 

Bubble 

Point 

Pressure 

1500 

1500 
1200 
900 

1500 

1500 

1500 
1200 
900 

900 

900 

900 

900 

900 

Porosity 

0.  35 

0.  35 

0.  25 

0.  25 

0.  25 

0.  25 

0.  25 

0.  20 
0.  25 
0.  30 

0.  25 

0.  25 

Connate 

Water 

Saturation 

0.45 

0.45 

0.  30 

0.  30 

0.  30 

0.  30 

0.  30 

0.  30 

0.  25 
0.  30 
0.  35 

0.  30 

(m) 
Gas  Cap  VoL 
Oil  Vol. 

0.1 

0.1 

0.  5 

0.1 
0.  5 
0.9 

0.1 

0.1 
0.5 
0.9 

0.1 

0.1 

0.1 

0.1 

M 

Viscosity 

Ratio 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

0.  8M 

1.  0M 
1.  5M 

1.0 

1.0 

1.0 

Case  No. 

1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

Figure 

No. 

4.0 

4.1 

5.0 

5.1 

6.0 
6.1 

7.0 
7.1 

8.0 
8.1 

9.0 
9.1 

10.0 
10.1 

11.0 
11.1 

12.0 
12.1 

none 

%      A  N/N 


23.5  26.0    25  22.5 
28.6    28.6    26.5    26.5    26.0    28.0    23.6      23.6    23.6    23.6 

28.6  30.1    21.0  25.0 
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produced  directly  from  the  gas  cap.     Since  the  computer  program  did  not 
account  for  such  gas  production,   an  estimate  of  the  production  from  the  gas 
cap  was  made  and  production  data  were  corrected  accordingly.     This  non- 
ideal  situation  illustrates  the  flexibility  of  the  program.     Water  influx  is  not 
included,   but  may  be  reflected  by  altering  the  relative  permeability  curve. 
Changing  the  input  data  requires  some  effort,  but  the  effort  is  slight  when 
compared  with  that  required  to  carry  out  manually  the  Muskat  Prediction. 

When  the  predicted  results  were  analyzed  and  a  set  of  parameters 
selected,   it  was  necessary  to  consider  whether  or  not  the  set  of  parameters 
chosen  for  Reservoir  B  was  unique. 

Case  I  illustrated  that  no  set  of  parameters,   employed  with  porosity 
of  0.  30,   0.  35  and  0.  40,   would  satisfy  observed  production  data.     Consequently, 
when  Case  II  variables  were  employed,   it  was  quickly  apparent  that  only  a 
porosity  of  0.  25  gave  actual  observed  production  in  barrels  of  stock  tank  oil. 
N  was  not  used.     After  the  porosity  was  established,   only  the  Kg/K0  vs.   SQ 
curve,   which  was  positioned  by  S      »   0.  30,   gave  a  prediction  which  super- 
imposed on  the  observed  production  data.     Since  only  this  single  combustion 
of  Sw  and  <t>  resulted  in  observed  production  (%     AN/N  vs.   Pressure)  which 
superimposed  on  predicted  production,  the  solution  is  a  unique  solution.   These 
two  variables  were  of  primary  importance  because  they  were  used  to  determine 
the  observed  %      A  N/N.     The  computed  %      A  N/N  is  a  function  of  (1  -  Sw)  times 
1/B0  and  does  not  reflect  porosity  at  all.     It  follows  that  when  an  Sw  satisfies 
both  the  observed  and  predicted  %     AN /N,   it  is  a  unique  value  of  Sw. 
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CONCLUSION 

In  this  paper,   the  IBM- 701  digital  computer  solution  of  the  Muskat 
equation  was  used  to  predict  a  solution  gas  drive  reservoir.     The  conclusions 
drawn  from  this  analysis  are  as  follows: 

1.  When  using  the  Muskat  Prediction  to  analyze  a  reservoir,   one 
prediction  is  usually  insufficient. 

2.  A  match  between  observed  and  predicted  data  is  obtained  only 
after  several  calculations.     During  these  calculations,   the 
different  parameters  were  varied  independently  and  in  groups . 

3.  Even  when  a  match  is  obtained,  the  effect  of  variation  of  the 
remaining  parameters  must  be.  investigated  to  be  reasonably 
sure  that  they  will  not  destroy  the  match. 

4.  The  Pressure  vs.  %     A  N/N  curve  and  the  GOR  vs.  %    A  N/N 
curve  could  not  be  simultaneously  matched  with  any  set  of 
reservoir  parameters.     The  discrepancy  between  the  observed 
and  predicted  GOR  curves  is  attributed  to  gas  production  directly 
from  the  gas  cap. 

5.  The  reservoir  engineer's  experience  and  knowledge  of  a  reservoir 
are  a  necessary  part  of  the  procedure.     He  must  decide  when  the 
parameters  best  specify  the  reservoir.     Since  only  the  ideal 
reservoir  will  conform  exactly  to  the  predicted  solution,   the 
reservoir  engineer  must  decide  when  a  reasonable  approximation 
has  been  made  and  stop  his  analysis  at  this  point. 
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NOMENCLATURE 


English 


Bo 
d 

f 

Gi 

Gp 

i 

K4 

m 

n 

Np 

P 

Pab 

Pb 
P 

Q 
R 

R, 


;g 

»o 

>oi* 


con 


U 
X 
Y 


oil  formation  volume  factor,   reservoir  barrels  per  stock  tank  barrel 

differential  operator 

fraction  of  produced  gas  injected  while  producing  below  bubble-point 

pressure  with  declining  pressure 

cumulative  gas  injected,   thousands  of  cubic  feet  of  gas  measured  at 

standard  temperature  and  pressure  per  acre-ft.   of  formation 

cumulative  gas  produced,   thousands  of  cubic  feet  of  gas  measured  at 

standard  temperature  and  pressure  per  acre-ft.   of  formation 

subscript,   indicates  initial  or  starting  condition  of  pressure  or 

saturation 

effective  permeability  to  gas,   md 

effective  permeability  to  oil,   md 

ratio  of  initial  reservoir  free  gas  volume  to  initial  reservoir  oil 

volume 

number  of  pressure  or  saturation  steps 

cumulative  oil  produced,  barrels  of  tank  oil  per  acre-ft.   of  formation 

pressure,   psia 

abandonment  pressure,    psia 

bubble-point  (saturation)  pressure,   psia 

constant  pressure,   psia 

see  Equation  3 

producing  GOR,   cubic  feet  of  gas  measured  at  standard  temperature 

and  pressure  per  stock  tank  barrel  of  oil 

solution  GOR  (gas  solubility  in  oil),    cubic  feet  of  gas  measured  at 

standard  temperature  and  pressure  per  stock  tank  barrel  of  oil 

gas  saturation,   fraction  of  pore  space 

oil  saturation,   fraction  of  pore  space 

oil  saturation  in  the  gas  cap,   fraction  of  pore  space 

water  saturation,   fraction  of  pore  space 

see  Equation  3 

see  Equation  3 

see  Equation  3 


Greek 


A   - 

ft' 


relative  gas  density,    cubic  feet  of  gas  measured  at  standard 

temperature  and  pressure  per  cubic  foot  of  gas  measured  at 

reservoir  temperature  and  pressure 

difference  (  P  x  PI  ~  Po  »   S0 

gas  viscosity,   cp 

oil  viscosity,   cp 

porosity,   fraction 


S2  "  si  > 


CALCULATION  OF  INITIAL  OIL  IN    PLACE 
Sand  count  of  Figure  3:     2520  acre  feet. 

Case  I : 

N    .       7758  (2520)  (<j>)  (1  -  Sw) 
Boi 

1.     B.p  =1500,   <t>  •   0.35,   Sw  »   0.45,   Boi    -  1.185 


N    ,(7758)  (2520)  (0.35)  (0.55)         .   ,.  „  ,  1Q6  „  bM, 

1.185 


2.     B.P*    900,  <t>  -   0.  35,   Sw  -   0.  45,   BQi  -  1. 113 


N    .(7758)  (2520)  (0.35)  (0.55)        ,   3.  38  x  106  ST  bbls 
1.113 


Case  II: 


3.     B.P^1500,   <b  '   0.25,   Sw  =   0.30,    Boi    »  1.185 


N    ,<7758)  (2520)  (0.25)  (0.7)       .   2.  88  x  106  ST  bbls 

1       lob 


4.     B.P-   900,   i    -   0.25,    Sw  =   0.30,    BQi  «  1.113 


M        (7758)  (2520)  (0.  25)  (0.7)         „,  6 

N    "   1  U3      — '  3.12x10    ST  bbls 


ii 


CALCULATION  OF  %  A/V//V 


Pressure                   bbjs     of                 Case  I       Case  II       GOR  GOR 

psia STO  Produced % % (actual)     (corrected) 


1100 

50,  000 

1.  48 

1.6 

1000 

100,000 

2.96 

3.  2 

900 

175,000 

5.18 

5.7 

800 

295,000 

8.  75 

9.5 

700 

410,  000 

12.1 

13.1 

600 

490,  000 

14.5 

15.7 

600  250 

900  500 

1200  800 

1670  1050 


Case  I:    N    ■  3.  38  million  bbls  of  STO,   <t>  *  .35,  Sw  =   .45 
Case  II:    N    »  3. 12  million  bbls.  of  STO,   4=  .25,  Sw  *   .30 


iii 
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